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Abstract: 

Hadronic events produced in e+e^ collisions by the LEP collider and recorded by the OPAL 
detector were used to form distributions based on the number of reconstructed jets. The data 
were collected between 1995 and 2000 and correspond to energies of 91 GeV, 130-136 GeV and 
161-209 GeV. The jet rates were determined using four different jet-finding algorithms (Cone, 
JADE, Durham and Cambridge). The differential two-jet rate and the average jet rate with 
the Durham and Cambridge algorithms were used to measure as in the LEP energy range 
by fitting an expression in which O^a^) calculations were matched to a NLLA prediction and 
fitted to the data. Combining the measurements at different centre-of-mass energies, the value 
of as{Mz) was determined to be 

as(Mz)=0.1177 ± 0.0006(stat.)±0.0012(expt.)±0.0010(had.)±0.0032(theo.) . 
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* Deceased 

1 Introduction 

In the Standard Model of elementary particle interactions, the strong interaction is described by the 
theory of Quantum Chromodynamics (QCD), and depends on just one fundamental parameter, the 
strong coupling ag. The value of as is expected to depend on the energy scale of the interaction. It 
is therefore an important test of the theory to determine the value of ag experimentally at as many 
different energies as possible. It is also important to use as many different techniques as possible, 
as different measurements of as are sensitive to different theoretical and hadronization variations. 

Indeed, many methods have already been employed to evaluate [1]. At very low energies the 
value of as can be measured using the hadronic decays of the r lepton and heavy quarkonia. Low 
energy determinations are also available using scaling violations and sum rules from deep inelastic 
scattering experiments. Higher energy determinations of ag come from collider experiments (e"'"e~, 
pp, pp or ep) using properties of the created hadron system which are explicitly dependent on the 
value of as{Q), where Q corresponds to the energy scale at which the interaction takes place^. 

During the LEP1.5 {^/s ~133 GeV) and LEP2 (above W+W" threshold) operational phases of 
the Large Electron-Positron collider at CERN, events were recorded with centre-of-mass collision 
energies ranging from 91 GeV to 209 GcV. Events of the form c+c^ ^liadrons can be used to deter- 
mine distributions based on the ensemble of final state hadrons (event shapes) or on the ensemble 
of jets (jet rates). Previous results by OPAL for an as determination based on event shapes and jet 
rates using the Z dataset collected during the LEPl phase can be found in [2]. Determinations of 
ttg from LEPl. 5 and LEP2 datasets up to 189 GeV have already been reported by OPAL based on 
event shape distributions [3, 4, 5] and on jet rates [6]. Another OPAL paper [7] uses the same data 
that have been presented here to measure event shapes. 

For the analysis presented in this paper we used data collected during the LEPl. 5 and LEP2 
phases to construct jet rate distributions using several jet clustering algorithms. The differential 
two-jet rate, D2, and the average jet rate, {N), were used to determine values of as{^/s) at the 
four combined centre-of-mass energies composed of data within the LEPl. 5 and LEP2 datasets. 
Theoretical predictions were fitted to these distributions to extract the value of as{^/s). 

The paper is organized as follows. Section 2 contains a brief description of the OPAL detector. 
A summary of the data and the Monte Carlo samples used in the analysis is given in Section 3. In 
Section 4, we define the jet rate distributions. The methods used to select signal events and reject 
backgrounds are presented in Section 5. The variations used for systematic studies are detailed in 

^For e"'"e~ collisions Q =-\/s. 
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Section 6. Finally, the results of this analysis are given in Section 7, followed by a conclusion and 
summary in Section 8. 

2 The OPAL Experiment 

A full description of the OPAL detector can be found in [8] . The critical components of the detector 
in the identification of jets were the central tracking chambers, which were used to reconstruct 
charged particles, and the electromagnetic calorimeters, which measured the total energy deposited 
by electrons and photons. 

The tracking chambers were located inside a solcnoidal magnet which provided a 0.435 T axial 
magnetic field along the beam axis. The main component of the tracking system was a large- 
volume jet chamber, which was approximately 4.0 m long with an outer radius of 1.85 m. The jet 
chamber was separated into 24 sectors, each with a radial plane of 159 sense wires separated by 1 
cm. The momenta of tracks in the x — y plane^ were measured with a precision parametrized by 
ap/p = v'0.022 + (0.0015 •p[GeV/c])2. 

The calorimetry systems were outside the solenoidal magnet. The electromagnetic calorimeter 
was composed of 11704 lead glass blocks in the barrel and endcap regions, representing about 25 
radiation lengths in the barrel and more than 22 in the endcap. The iron sampling hadron calorimeter 
was located just outside the electromagnetic calorimeter, and provided the stopping power to contain 
most hadronic showers. Luminosity was determined using small-angle Bhabha events detected in 
the forward detectors and silicon-tungsten calorimeter [9]. 

After an event was triggered [10], data were collected from the subdetectorsand processed by 
the OPAL data acquisition system [11]. The raw event data were transferred to a farm of computer 
processors where the events were fully reconstructed and written to tape for offline analysis. 

3 Data and Monte Carlo Samples 

The data used in this analysis were collected by OPAL between 1995 and 2000 and correspond to 
integrated luminosities of 14.7 pb~^ of data taken with centre-of-mass energy 91 GeV, 11.3 pb~^ of 
LEP1.5 data with centre-of-mass energies between 130 GcV and 136 GeV and 707.4 pb~^ of LEP2 
data with centre-of-mass energies ranging from 161 to 209 GcV. The 91 GeV data, known as Z- 
calibration data, were primarily collected for calibrating parameters used in the OPAL reconstruction 
algorithms. This Mz sample had the same detector configuration as the other centre-of-mass energy 
points. The exact breakdown of the centre-of-mass energies together with the respective luminosities 
and numbers of selected events are given in Table 1. The thirteen points in Table 1 represent the 
main samples of the spread of energies in the LEP1.5 and LEP2 data. 

The data were combined into four datasets. The LEP1.5 data provided a single energy point 
at an event-weighted centre-of-mass energy of 133 GeV, while the LEP2 data were split into two 
energy points, one with an event-weighted centre-of-mass energy of 177 GeV using data in the range 
161-185 GeV (with a total integrated luminosity of 78.1 pb~^) and another at 197 GeV (with a 
total integrated luminosity of 628.3 pb~^) using data in the range 188-209 GeV. Together with the 
Z-calibration data this provided for a determination of ccg at four centre-of-mass energies. 

^The right-handed OPAL coordinate system is defined so that z is the coordinate parallel to the e~ beam direction 
and the x axis points to the centre of the LEP ring, r is the distance normal to the z axis, is the polar angle with 
respect to the z axis and <p is the azimuthal angle with respect to the x axis. 
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A number of Monte Carlo samples were created to correct for detector acceptance and resolution 
effects, to correct for hadronization effects and to estimate the contribution of background processes. 
These Monte Carlo samples were produced using a full simulation of the detector [12], followed by 
the same reconstruction and selection algorithms applied to the real data, and are referred to as 
"detector-level" samples. Other samples without the full detector simulation are discussed in Section 
5.2. 

PYTHIA 6.150 [13] was used to provide the default Monte Carlo samples (for the process 

e'''e~ Z/7* ^ qq — ^ hadrons) which were used to correct the high energy datascts. The Z- 
calibration dataset was corrected using JETSET 7.408 [14]. The use of JETSET for the lower 
energy data is a matter of convenience only, and not due to any inconsistencies in PYTHIA at 
this energy. Any differences between the two generators is expected to be negligible. Hadroniza- 
tion corrections were evaluated by comparing results with an alternative Monte Carlo sample, 
HERWIG 6.2 [15] which uses the cluster model of hadronization. This was compared with the 
string model of hadronization in PYTHIA. The parameters which were involved in the Monte Carlo 
simulation, both for JETSET/PYTHIA and HERWIG, were tuned to OPAL data collected at the 
Z peak, including global event shapes, particle multiplicities and fragmentation functions [16, 17]. 
The generation of the initial quark-antiquark pair for each Monte Carlo sample was implemented at 
LEP2 using the /C/C2f 4.13 event generator [18], which has an improved description of photon pro- 
duction in the initial and final states with respect to the one currently implemented in the PYTHIA 
generator. The available detector-lcvcl Monte Carlo samples are listed in Table 1. 

Above the W^W^ production threshold (161 GeV), the main background was expected to come 
from four-fermion events (e+e~ ^W^W~^4f), in particular those events in which two or all four 
of the fermions were quarks. The contribution of these backgrounds in data was estimated using 
Monte Carlo samples generated using KORALW 1.42 [19] (for qqq'q' and qqil^'^ where £ = e, ij,,t,u 
but a 7^ e+e^) and grc4f 2.1 [20] (for eeqq). Grc4f 2.1 was used to generate all the expected 
four-fermion background samples for the 161 and 172 GeV data. The background distributions 
were normalized to the luminosity of the dataset and subsequently subtracted from the measured 
distributions. The LEP1.5 energies were well below the W'''W^ and ZZ production thresholds [21] 
and were therefore expected to have no significant four-fermion backgrounds. The total expected 
background contribution from "four-fermion" e^e^-^ qqqq events is 1.2% of the combined LEP1.5 
data sample and it was neglected in the analysis. 



4 Jet Rate Distributions 

Jets were formed from the final state objects by applying jet clustering algorithms. These algorithms 
use the kinematic and spatial (geometric) properties of the individual objects in order to classify 
them as belonging to a specific jet. We used here the Durham [22], Cambridge [23], JADE [24] and 
the R and e variants of the Cone [25] jet clustering algorithms. 

The Durham and Cambridge algorithms construct a test variable built from the energy and 
angular separation between two particles, 

_ 2uim{Ef,E]}{l - cos%) 

vis 

where Ei is the energy of particle i, 9ij the angle between the particle i and j and E^^g is the total 
visible energy in the event. The pair that produces the smallest value of yij is chosen first. The value 
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of this test variable is compared to a predefined parameter, z/cut) called the jet resolution parameter. 
If the test variable is smaller than j/cut particles i and j are merged into a pseudo-particle. Merging 
means that the momenta of particles i and j are removed from the set of momenta and the the sum 
of their four-momenta is added to the set of momenta. After the merging, the clustering starts again 
using the momentum set and it continues until all test variables become larger than ?/cut- After the 
clustering stops, all remaining (pseudo-) particles are classified as jets. 

The Cambridge algorithm differs slightly from Durham in its implementation. In the Cambridge 
algorithm particles are first paired together by minimizing the variable Vij = 2(1 — cos dij). The 
standard test variable is then constructed and compared to the jet resolution parameter, ycut- The 
procedure followed is then identical to that of the Durham algorithm, except that Cambridge freezes 
out soft jets by accepting only the lowest energy (pseudo-)particle as the jet when yij > j/cut- The 
number of jets reconstructed in the event, using Diirham or Cambridge, is therefore a function 
of the jet resolution parameter. The JADE algorithm follows the same procedure as the Durham 
algorithm; however, it uses the scaled invariant mass, yij = 2EiEj(l — cos 6ij)/E^^^, of particles i 
and j as the test variable. 

In the Cone jet finding algorithm, a jet is defined as a set of particles whose three-momentum 
vectors lie inside a cone of half angle R, where the direction of the sum of their three-momentum 
vectors defines the cone axis. In addition, the total energy of the particles assigned to a jet is 
required to exceed some minimum value e. Typical values are R = 0.7 rad and e = 7 GeV for jets 
in e"'"e~ annihilation at LEP 1 energies. When analysing events at the detector level, we replaced e 
hj e' = e ■ E^is/y/s to compensate for the incomplete detection of the energy of the event. In our 
studies, the jet rate was computed at fixed s = 7 GeV as R was varied, and at fixed R = 0.7 as e 
was varied. The former is sensitive to the angular structure of jets, and the latter to their energy 
distribution. 

The fraction of multihadronic events in a given sample that are classified as containing n jets 
for a given value of the jet resolution parameter (ycut; R or e) is referred to as the ra-jet rate. This 
n-jet rate is explicitly defined as 

Mvcu.) = - (1) 

ftot ^Vtot 

where a"„ is the cross-section for the production of a hadronic event with n jets at fixed ycut, <^tot is 
the total hadronic cross-section, Nniycnt) is the number of events in a sample with n jets for a given 
value of ycut and Ntot is the total number of events in that sample. 

The differential n-jet rate was also determined. It is the derivative of the n-jet rate with respect 

to J/cut, 

DM = (2) 

For the case when n = 2, the differential 2-jct rate reduces to D2 = y23, where 7/23 is the value of 
the jet resolution parameter where the event flips from a 2- to a 3-jet event. When the Durham 
algorithm is used to define jets the value of D2 (denoted y^) is also an event shape variable. 

The average number of jets per event in a given sample, as a function of the jet resolution 
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parameter, is defined to be 

imiycnt) 



A QCD prediction which matches an ©(ctg) (next-leading order) prediction, based on the QCD 
matrix elements [26], with a resummed, next-leading logarithmic approximation (NLLA) [27] pre- 
diction, such that terms that appear in both predictions are not double counted, was fitted to data. 
In this analysis we used the InR matched D2 [28, 29, 22] and (N) [6, 30] predictions to fit to the 
distributions of the observablcs. The difi"crcntial and average jet rates were determined using the 
Durham and Cambridge clustering algorithms, since resummed predictions only exist for these al- 
gorithms. This provided four separate observables {D2, D2 , {N)^ and (iV)*^) which were used to 
determine a value of Og at the four different centre-of-mass energy values. 

5 Analysis Procedure 

5.1 Selection Method 
5.1.1 Preselection 

All events within a dataset were required to contain information from both the central jet chamber 
and the electromagnetic calorimeter, meaning both these subdetectors must have been flagged as 
being on and in good operational condition. In addition events were required to be tagged as 
multihadronic in order to be preselected for analysis. Multihadronic events were identified using 
the criteria described in [31] for events with >Mz and in [32] for Z-calibration events. To pass 
the preselection, an event was required to contain at least seven good tracks to reduce potential 
backgrounds arising from the production of r leptons (e+e~ — r+r") decaying into hadrons and 
from two-photon interactions producing quarks. Good tracks were defined as those which had 

• at least 40 hits in the jet chamber 

• at least 150 MeV/c transverse momentum relative to the beam axis. 

• the distance of closest approach to the interaction point in the r — cf) plane satisfying do < 2 cm 

• the point of closest approach < 25 cm from the interaction point in the z-direction 

Clusters of energy in the calorimeters were also used in the analysis; good clusters were defined as 
those which produced a signal in at least one block in the barrel electromagnetic calorimeter corre- 
sponding to an uncorrected energy of at least 100 MeV or of 2 blocks in the endcap electromagnetic 
calorimeter corresponding to an uncorrected energy of 250 MeV. The hadron calorimeter was not 
used in this analysis. 

All of the good quality tracks and clusters in the event were used to define "objects" representing 
particles using an algorithm (MT) to correct for double counting of energy. This MT algorithm 
produced a uniquely defined array of track and cluster objects. The trajectories of the tracks 
measured in the central tracking chambers were extrapolated to the clusters in the electromagnetic 



— V'nc7„(ycut) 

— ^nNnivcMt)- 



(3) 
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calorimeters. If the energy of the cluster was less than expected from the track, then the cluster was 
omitted to avoid double counting of energy, since the momentum resolution for tracks was typically 
better than the calorimeter energy resolution. If the energy of the cluster was larger than expected 
the energy of the cluster was reduced by the expected amount with the remaining energy interpreted 
as due to photons or neutral hadrons. These remaining clusters and those which were not matched 
defined the four-vectors of "neutral" particles. In all cases tracks were treated as charged pions and 
neutral particles were treated as being massless. 

5.1.2 Containment 

We ensured that most particles in the event were well contained in the detector and not lost down 
the beam line by imposing a cut on the direction of the thrust axis [33] , 

• IcoseiTl < 0.9, 

where is the angle between the beam axis and the direction of the thrust axis. The thrust axis 
direction was determined from all tracks and clusters in the event, without correcting for double 
counting with the MT algorithm. 

5.1.3 Initial State Radiation (ISR) Cuts 

The events of interest for this analysis were e+e""— qq events where the final-state qq pair had 
the full centre-of-mass energy. The effective centre-of-mass energy of the e+e" collision can be 
reduced by the emission of one or more ISR photons. At LEP2, approximately three quarters of the 
multihadronic events were such "radiative return events" , where the invariant mass of the qq pair 
was close to the Z mass. The effective centre-of-mass energy of the collision after ISR, a/s' [34], was 
evaluated, and the requirement 

• y^-Vs' <10 GeV 

was imposed to select full-energy events. 

To calculate s/s', all isolated photon candidates with energies greater than 10 GeV were iden- 
tified. The Durham jet reconstruction algorithm [31] was then used to group the remaining tracks 
and clusters into jets. ISR photons arc often emitted close to the beam direction. Three kinematic 
fits were performed, under the assumptions that 

• there were two undetected photons (in opposite directions along the beam pipe) 

• there was one undetected photon 

• all photons were observed in the detector, 

respectively. The fit with the most acceptable was selected, and s/s' was calculated from the 
invariant mass of the jets, excluding any photons. 

The power of this cut can be seen in Figure 1. The efficiency for selecting non-radiative qq 
events is given in Table 2. The purity of non-radiative events was found to be approximately 73% 
in all of the LEP1.5 and LEP2 data samples. Non-radiative qq events are defined as those in which 

— ^/s[^ < 1 GeV, where s^rue determined from generator-level information in the PYTHIA 
samples. This ISR cut was applied to all analyzed datasets with the exception of the Z calibration 
data. 
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5.1.4 Final Cuts 



The dominant background to the process e"'"e~ ^ Z/7* ^ qq ^ hadrons at LEP2 came from the 
four-fermion process e^e~ ^W^W~ in which one or both of the bosons decayed hadronicahy, pro- 
ducing two or four quarks in the final state. This background was expected to make up approximately 
30% of all observed events which pass the first stage of cuts in each of the LEP2 datasets. These 
backgrounds were addressed by placing a cut [35] on two likelihood values which indicate how likely 
an event is to be a non-QCD four-quark or a semi-leptonic event: 



• C,^lv < 0.50 

The effect of these cuts in each of the LEP2 datasets and the expected backgrounds can be seen on 
Figure 2. 

The four-quark likelihood value [36], ^qqqq, was estimated from four kinematic variables de- 
scribing characteristics of hadronic W+W~ decays like their four-jet nature and angular structure. 
These variables were used to construct event probabilities based on two hypotheses: first, that the 
event was due to a hadronically decaying W"'"W~ pair (W"'"W~^qqqq) and, second, that the event 
was due to a hadronically decaying Z/7* (e"'"e~ Z/7* ^qq). The probabilities were combined 
to produce the discriminating likelihood, -Cqqqq- This cut reduced the expected background by 
approximately 80% so that it constituted only 9% of the observed number of events. 

The semi-leptonic likelihood [36], Cqqi^, was based upon three separate likelihoods, one for each 
lepton species (i =e, /i, r). Each of these likelihoods was based on ten variables describing the 
properties of the lepton, the jets produced by the qq pair and the missing energy carried away by 
the neutrino. This cut in conjunction with the cut on the four-quark likelihood removed almost 90% 
of the background expected in the observed LEP2 dataset. The effect of these cuts can be seen in 
Figure 2. 

These likelihood cuts also reduced the backgrounds arising from e+e" ^ZZ in which one or 
both of the Z bosons decay hadronically. ZZ production contributed only a small fraction of the 
background due to its lower cross-section compared to W+W~ production in the energy ranges used 
in this experiment. The likelihood cuts were applied only to those datasets with >161 GeV. 

The expected size of the total background contribution to each dataset was determined by Monte 
Carlo predictions after scaling to the luminosity of the dataset. The effect of the final cuts and the 
expected four-fermion backgrounds for each centre-of-mass energy dataset can be seen in Table 2. 
As seen in the table the likelihood cuts greatly increased the purity of selected non-radiative qq 
events. The LEP2 datasets data with ^/s > 183 GeV were typically ~70% pure following the ISR 
cuts; however, after the final cuts this increased to a purity of 94-95%. 

5.2 Monte Ccirlo Corrections 

The values of the variables Rn, Dn and (iV) were determined for each accepted event using the 
MT-corrected tracks and clusters. These values were then compiled into histograms as a function 
of VcvA with bins of varying size. The background rejection cut did not completely remove all of the 
expected background events from W''"W~ and ZZ production (referred to as four-fermion background 
in this paper). The remaining backgrounds, taken from the Monte Carlo, were subtracted from 
the corresponding measured distributions on a bin-by-bin basis. Systematic uncertainties in this 
procedure will be discussed in Section 6. 
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Corrections to the distributions were also made for effects arising from finite detector resolution 
and a limited detector acceptance (recall that the cut on | cos 6't | reduced the fiducial volume) 
and for residual ISR events which were not removed by the \fs' cut. These corrections were done 
separately for each variable and were accomplished by comparing distributions from two separate 
Monte Carlo samples, one of which had gone through a full detector simulation, including effects 
of detector resolution and acceptance and initial state radiation, called the "detector" level. The 
other sample used only the generator-level hadrons and had not gone through the detailed detector 
simulation. In this sample all short-lived particles (r < 3 x 10~^''s) had decayed and a requirement 
that \fs — \/s[^^ < 1 GeV was imposed. This "hadron level" sample was thus expected to produce 
distributions arising solely from the properties of the underlying hadrons, free of any detector biases 
determined over the full acceptance without any limitations arising from limited resolution. 

Correction factors for each bin of the distributions were determined from the ratio of the two 
Monte Carlo distributions. Thus, any bin, i, of the measured distributions was corrected via 

/-i/MC\ 

^f*' = J (^'^'^ - (4) 

where 7i and V represent distributions at the hadron level and the detector level respectively, and 
D^^'^^'^corresponds to the expected size of the total background in bin i. The hadron level was used 
in this analysis when determining jet rate distributions; theoretical predictions which wore fitted 
to these distributions were obtained from computations valid at the "parton level". The parton 
level corresponds to distributions that would be produced if only the partons created immediately 
following the e"'"e~ annihilation and before the hadronization phase were used in the analysis. The 
parton level Monte Carlo sample was biiilt from quarks and gluons that were produced during the 
parton shower simulated by the generator before the hadronization phase began. As in the case of 
the hadron level, the parton level sample gave rise to distributions that were free of initial state 
radiation without any detector simulation applied. The correction factor determined from the ratio 
of the hadron level to the parton level was applied to the theoretical predictions before the fitting 
procedure. This factor corrected the prediction to the hadron level so that it could be compared to 
the corrected hadron level distribution determined from the data, 

where H and V represent distributions at the hadron level and the parton level respectively. Sta- 
tistical errors on the data distributions included the effects of Monte Carlo statistics. 

6 Systematic Uncertainties 

6.1 Experimental Systematic Variations 

Several selection algorithms and selection cuts were varied to determine their impact on the results 
of the analysis. In all cases the result from the variation was compared to the result from the 
standard selection, and the difference was taken as a contribution to the total systematic error. The 
systematic variations that were used (given in descending order of their contributions to the overall 
experimental uncertainty) are: 
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Track-Cluster Matching Systematic errors relevant to the definition of objects defined witliin 
the tracking chambers and calorimeters, and hence used for jet clustering, were estimated by 
comparing the results using the MT package to a method using all selected tracks and clusters 
without taking into account the possibility of double counting. 

Detector Correction The uncertainty in modelling the detector was investigated by using HERWIG 

Monte Carlo datasets in place of PYTHIA to correct from the detector level to the hadron level. 
The hadronization correction was still performed using the PYTHIA Monte Carlo samples. 

Containment The constraint on the direction of the thrust axis was tightened to |cos^t| < 0.7, 
restricting events to the barrel region. 

ISR A possible systematic effect introduced through the selection of events with little or no ini- 
tial state radiation was estimated by repeating the selection using a second method of ISR 
determination [37]. This alternative method assumed while performing the kinematic fit that 
there was always a single photon which either escaped undetected down the beam line or was 
detected in the electromagnetic calorimeter. 

Cqqqq and jCqqii, To account for the systematic uncertainty that arises from the value of the cut 
placed on the W"'"W~ hadronic and semi-leptonic likelihoods, the values of the likelihood cuts 
were changed to 0.1 and 0.4 for Cqqqq and to 0.25 and 0.75 for Cqqii,. In each case the largest 
deviation was taken as the systematic uncertainty. 

Backgrounds To account for uncertainties introduced during background subtraction, arising from 
imprecise knowledge of the four-fermion cross-sections, these cross-sections were conservatively 
varied by ±5% and the largest deviation from the standard value was used to determine an 
overall systematic error. 

It should be noted that there was no single dominant contribution to the overall experimental 
systematic error. In general, the largest contributions occurred for the track-cluster matching and 
detector correction variations, while the least significant contributions came from the variations of 
the background cross-section and the semi-leptonic likelihood. 

6.2 Hadronization Systematic Variations 

Systematic uncertainties arise from the modelling of the hadronization. These were estimated by 
using HERWIG Monte Carlo samples, which employ a different hadronization model, in the deter- 
mination of the theoretical prediction correction factor (Eqn 5). The PYTHIA-HERWIG differences 
were taken to be the systematic uncertainties from hadronization. The statistical component due 
to limited Monte Carlo statistics was included in the total determination of the hadronization un- 
certianty. 

6.3 Theoretical Systematic Veiriations 

Three further systematic variations were considered when fitting the predictions to the data to 
determine as- 

Fit Range We investigated the choice of the range of bins used in performing the fit of the the- 
oretical prediction to the data. The fit range was first increased by two bins, by adding one 
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bin to each endpoint of the fit range. In the case where one of the endpoints was aheady 
at the maximum ahowed bin value, only the other point was extended by one bin. A second 
variation decreased the fit range by two bins, by removing a bin from each endpoint of the fit 
range. The largest deviation from the standard fit value was taken as the contribution to the 
systematic error. 

Renormalization Scale The second fit-related systematic variation accounted for the uncertainty 
due to the dependence of the fixed order and resummed predictions on the renormalization 
scale, x^, where = The value of X/^, which was set to 1.0 in the standard fits, was 

varied to 0.5 and 2.0 respectively. The largest deviation from the standard fit value was taken 
as the contribution to the systematic error representing renormalization scale uncertainty. 

Logarithm rescaling, a;_L In the resummation process for event shape variables, like D2, there is 
an arbitrariness due to the definition of the logarithms which were resummed. In this analysis 
we used ln(l/|/cut); however, this could be generalized to powers of ln[l/(xL ycut)] [38]. 
The standard value for this rescaling (xz,=l) was varied to | and |, to investigate the system- 
atic effect of this arbitrariness. No analoguous rescaling prescription has yet been developed 
for the case of the average jet rates, so the xl variations are only shown for the D2 distributions 
for comparison to the x^ variation. Hence this variation was not used in the determination 
of the total theoretical systematic error and was included only as a cross-check and for com- 
parison with other analyses. A comprehensive study of the combination of various theoretical 
variations to produce a global theoretical uncertainty for event shape distributions is given 
in [39]. 

The differences between the standard result and those due to the above variations were separated 
into three categories: experimental, hadronization and theory (see Table 3). The variations in each 
of these categories were added in quadrature and the result taken as the systematic uncertainty 
for that category. In the case of asymmetric errors, the error was symmetrized by taking the 
largest systematic variation and applying it as the full systematic contribution. A summary of the 
systematic variations is provided in Table 3. 

7 Results 

Data from thirteen datasets were used in this analysis: one Z-calibration dataset, two LEP1.5 
datasets and ten datasets from LEP2. These thirteen datasets were combined to produce four 
higher statistics datasets (-y/s =Mz, 133, 177 and 197 GeV) which were analysed separately. The 
raw distributions (n-jet fractions, D2 and (N)) for each of the datasets were determined as functions 
of the jet resolution parameters defined using four different jet clustering algorithms. The distribu- 
tions underwent a bin-by-bin correction which included subtraction of expected backgrounds and 
correction for detector and residual ISR eflFects. Systematic eflFects were examined by varying the 
parameters used in selecting events (see Table 3). The difference between the corrected distribu- 
tions using these variations and those from the standard selection then determined the size of the 
systematic error on each bin of the distribution. 

Matched predictions of NLO and resummed calculations were fitted to the corrected Durham and 
Cambridge D2 and (N) distributions over a predetermined fit range (see Section 7.2.1), taking into 
account bin-to-bin correlations in (N) (see Section 7.2.2). These fits provided four values of as with 
statistical and systematic errors at each of the four centre-of-mass energies. Taking into account the 
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statistical and systematic correlations between the four measurements, they were combined into a 
single as result at each energy. 

7.1 n-Jet Fractions 

The n-jet fractions for the Z-calibration sample and those for the LEP1.5 and the two LEP2 samples 
arc shown in Figures 3 to 7. Each plot shows the fraction of events in a given sample that were 
determined to be n-jets for a given value of the jet resolution parameter at the hadron level^. The 
jet fractions were calculated using four different clustering algorithms. For the Cone algorithm, 
results from using both the R and e resolution parameters are plotted, showing the individual n- 
jet fractions for n < 2, n = 3 and n > 4 in Figures 3 to 4. Similarly the n-jet fractions for the 
JADE algorithm are shown in Figure 5 for n = 2,3, 4, 5. Figures 6 to 7 show the individual n- 
jet fractions for n = 2, 3, 4, 5 for the Durham and Cambridge algorithms respectively. The error 
bars on the points represent the total statistical and systematic errors added in quadrature. The 
Monte Carlo expectations corresponding to PYTHIA and HERWIG are also displayed on each plot, 
for each algorithm and energy. The Monte Carlo expectations match the measured n-jet fractions 
reasonably well. 

The differential two-, three- and four-jet rates, Dn, are plotted as a function of j/cut for the 
Durham and Cambridge algorithms for the Z-calibration, the LEP1.5 and the two LEP2 data samples 
in Figures 8 and 9 respectively. Similarly, the average jet rates are plotted as a function of i/cut for 
the Durham and Cambridge algorithms for the Z-calibration, the LEP1.5 and the two LEP2 data 
samples in Figures 10 and 11 respectively. The curves on all the plots represent the expected Monte 
Carlo distributions. There is good agreement between the data and the expectations from both 
PYTHIA and HERWIG. 

7.2 Fits to Determine the Value of cXg 
7.2.1 Differential Two-jet Rates 

The range over which the Z?2 and {N) distributions were fitted was determined by splitting the 
91 GeV (and 189 GeV) PYTHIA Monte Carlo samples into 100 statistically independent subsamples. 
The QCD predictions for each distribution were fitted to each of these subsamples, with ag as a 
free parameter, for a number of possible end-points of the fit range (requiring that the fit range be 
at least six bins). A P^r degree of freedom was determined for each fit range. The values 
were averaged over the 100 subsamples. The fit range that produced the smallest average P^r 
degree of freedom was chosen to be the default fit range for the Z-calibration (and high energy) 
datasets. Where no clear x'^ minimum was found, the largest reasonable range was chosen. The 
size of the fit range was then adjusted to ensure the range did not extend into a region where the 
hadronization corrections exceeded 10%, in particular at smaller ?/cut values. A further adjustment 
was made to exclude fit ranges where one of the endpoints produced a contribution of more than 30% 
to the overall x^ value. Potential correlations introduced in the correction of the distributions were 
accounted by including a covariance matrix in the x^ fit. The covariance matrices were determined 
in the manner detailed in Section 7.2.2. 

The fits of the Ini? matching prediction to the D2 distribution are shown in Figure 12 for the 
Cambridge algorithm and in Figure 13 for the Durham algorithm. The values of as determined from 

■^Further details of the data will be made available in the HEPDATA database, 
http://durpdg.dur.ac.uk/HEPDATA. 
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the In R fits for tlie four datasets, together with a complete breakdown of statistical and systematic 
uncertainties, are given in Tables 4 and 5 for the Cambridge and Durham'^ algorithms respectively. 

7.2.2 Average Jet Rates 

For the average jet rates, (N), the statistical errors were strongly correlated between points, since 
the same events were used to determine (N) at each value of ycut- The fit was performed using a 
correlated fit in which the covariance matrix was determined from the PYTHIA Monte Carlo 
sample divided into many detector level subsamples. Each subsample was then corrected to the 
hadron level using a second, statistically independent, PYTHIA sample. There were 1500 subsamples 
created for the Z-calibration dataset and 1000 subsamples created for the high energy datasets. The 
number of subsamples was chosen so that the elements of the covariance matrix would be stable 
and would have negligible fluctuations. These subsamples were used to build a standard covariance 
matrix, which was then converted into a correlation matrix. The statistical errors on the bins of 
the data distributions were then applied to this correlation matrix to produce the covariance matrix 
used in the fits. 

An example of the fit of the In R matching prediction to the {N) distribution is seen in Figure 14 
for the Cambridge algorithm and in Figure 15 for the Durham algorithm. The values of ccg deter- 
mined from the In R fits for the four datasets, together with a complete breakdown of statistical and 
systematic uncertainties, are given in Tables 6 and 7 for the Cambridge and Durham algorithms 
respectively. 

7.2.3 Running of ccg 

The four values of ag determined from the D2 and {N) distributions were combined into a single 
value at each centre-of-mass energy. The large statistical correlations between the four values were 
handled in a manner similar to that for the bin-to-bin correlations of the average jet rates. For each of 
the four distributions, 1000 Monte Carlo samples were used to determine the statistical correlations 
between the CKg values. The correlation matrices determined for the four separate energy points are 
given in Table 8. Using the statistical error for the ag value from each observable, the statistical 
covariance matrix was then determined. The full covariance matrix also included contributions from 
experimental, hadronization and theoretical uncertainties 

V = Vstat + V^xpt + Vhadr + Vtheo- (6) 

A weight was determined for each of the ccg values from the inverse of the covariance matrix, 
Wi = ^j{V~^)ij/^-{V~^)ij. The combined ctg value was then determined from the weighted 
sum of the as values. Only statistical and experimental systematic uncertainties were allowed to 
contribute to the off diagonal elements of the covariance matrix V, to ensure undesirable features 
such as negative weights were avoided. Hadronization and theoretical systematics were added only 
to (lie diagonal ckMncMits of IIk^ covariance matrix. 

*Notc that [7] also determines Qs using , denoted there as y^j. The small differences between the results have 
been investigated in detail, and are not significant. They may be attributed to differences in fit regions, the use 
of statistically different Monte Carlo samples, and the adoption of slightly different strategies for the assessment of 
theoretical errors. 



14 



The statistical and experimental uncertainties on the combined ag value were determined from 
the product of the weights with the individual covariance matrices [40], 

''"err = w^VerrW where err = stat, expt. (7) 

The hadronization and theoretical systematic uncertainties were determined by repeating the com- 
bination for each systematic variation separately using the same weights. The difference between 
these systematic combinations and the central value is taken as the systematic contribution to the 
error on the central value. For the experimental systematic covariance matrix, off diagonal elements 
were determined using a "minimum overlap" method, 

{Ve^pt)ij = max[(yexpt)ii, {Ve^pt)jj]- (8) 

The values of as determined for each centre-of-mass energy are given in Table 9 along with the 
breakdown of the uncertainties, both statistical and systematic. A comparison of the combined ag 
values in Table 9 with those determined from the individual D2 and (A^) distributions is seen in 
Figure 16. Taking the combined ag values at each centre-of-mass energy as an input, the value of 
Og can be run back to the Z pole using an C(as) prediction. These as(Mz) values are also shown 
in Tabic 9 and plotted against the world average value of as(Mz)=0.1187±0.0020 [41] in Fig 17. 
Using these values and taking into account proper statistical and systematic correlations, a weighted 
mean of as(Mz)=0.1177± 0.0006(stat.)±0.0012(expt.)±0.0010(had.)±0.0032(theo.) is determined. 
The four combined as values are plotted in Figure 18 as a function of the centre-of-mass energy, 
compared to the 0(ag) energy evolution of ag based on the determined value of ag(Mz). 

8 Summary and Conclusion 

Data from twelve LEP 1.5 and LEP 2 datascts, with ccntrc-of-mass energies ranging from 130 GeV 
to 209 GeV, were combined into three higher statistics datasets. These datasets and one combining 
several Z-calibration runs at 91 GeV were used to determine the n-jet fractions, the differential n-jet 
rates and the average jet rates for each of the energies. The different jet multiplicity distributions 
were compared to both PYTHIA and HERWIG Monte Garlo expectations. 

Hadron-lcvcl n-jet fractions were determined using four jet-clustering algorithms. Cone, JADE, 
Durham and Cambridge. For the Cone algorithm, measurements of the fraction of events with n < 2, 
3, >4 jets were presented as functions of R and e. In the case of JADE, Durham and Cambridge, 
measurements of the 2-, 3-, 4-, and 5-jet fractions were presented as functions of the jet resolution 
parameter, j/cut- In all cases there was generally good agreement between the measured jet fractions 
and the Monte Carlo expectations. 

Hadron-level determinations of the differential 77,-jet and average jet rates were performed for 
the Durham and Cambridge algorithms. The differential two-jet rate, D2, and the average jet rate, 
(N) were used to determine the value of as(\/s) for each of the four combined energy points. The 
determinations were carried out by fitting the InR matching predictions, appropriately corrected 
to the hadron level, to the hadron level data distributions over an appropriate range of j/cut; with 
as{\fs) as the variable parameter. The running of as(-y/i) was demonstrated by comparing the four 
values of ag as determined from the combined datasets as a function of their centre-of-mass energies. 

Using the measured values of as{^/s) a value of ag(Mz) was determined for each of the four 
datasets. A weighted mean taking account of correlations determined a final value of the strong 
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coupling at ^/s =Mx of 

as(Mz) = 0.1177 ± 0.0006(stat.) ± 0.0012(expt.) ± 0.0010(had.) ± 0.0032(theo.) 

where the error contains contributions from statistical, experimental, hadronization and theoretical 
uncertainties. The error on the determined value is slightly larger than that for the previous OPAL 
publication [6] which also used resummed predictions for D2 and average jet rate distributions, 
but explored slightly different energy ranges, including 35 and 45 GeV, and all LEP energies up 
to only 189 GeV. There is good agreement between these four values of as{Mz) measured in this 
analysis and previous determinations of as summarized in [1] and with the world average value of 
0.1187±0.0020 [41]. 
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Energy (GcV) 


Integrated 


Number of Events 


Year 


Nominal 


Range 


Mean 


Luminosity 

(pb-i) ' 


Data 


JETSET/ 

PYTIIIA 


HERWIG 


1996-2000 


91 


91.0-91.5 


91.3 


14.7 


426194 


459k 


443k 


1995, 1997 


130 


130.0-130.3 


130.1 


5.3 


1628 


50k 


50k 


1995, 1997 


136 


135.7-136.2 


136.1 


6.0 


1504 


50k 


50k 


1996 


161 


161.2-161.6 


161.3 


10.1 


1369 


100k 


80k 


1996 


172 


170.2-172.5 


172.1 


10.4 


1285 


100k 


80k 


1997 


183 


180.8-184.2 


182.7 


57.7 


6027 


100k 


100k 


1998 


189 


188.3-189.1 


188.6 


185.2 


18559 


500k 


100k 


1999 


192 


191.4-192.1 


191.6 


29.5 


2866 


200k 


100k 


1999 


196 


195.4-196.1 


195.5 


76.7 


7076 


200k 


100k 


1999, 2000 


200 


199.1-200.2 


199.5 


79.3 


6676 


200k 


100k 


1999, 2000 


202 


201.3-202.1 


201.6 


37.8 


3225 


200k 


100k 


2000 


205 


202.5-205.5 


204.9 


82.0 


6721 


200k 


100k 


2000 


207 


205.5-208.9 


206.6 


138.8 


10987 


375k 


100k 



Table 1: Integrated luminosity and the total number of preselected events for all samples. Included 
in this table are the Monte Carlo samples used to correct and make comparisons withthe data. The 
91 GeV dataset corresponds to the data collected during the Z-calibration runs. 
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Dataset 


ISR Fit 


Final 


Dataset 


ISR Fit 


Final 


91 GcV 




395695 


395695 


192 GeV 


Data 


717 


514 


MC Expected 


696±14 


471±4 


130 GeV 


Data 

H IT 1 


318 


318 


non-rad eff(%) 


82.8±0.7 


77.4±0.7 


MC Expected 


368 ± 3 


368 ± 3 


4f bkg frac (%) 


31.3±1.5 


5.3±1.2 


non-rad en(7oj 


85.4 ± 1.4 


85.4 ± 1.4 


196 GeV 


Data 


1732 


1137 


136 GeV 


Data 




'^^ 9 

OIZ 


MC Expected 


1746±24 


1162±9 


MC Expected 


329 ± 3 


329 ± 3 


non-rad eff(%) 


82.7±0.7 


77.1±0.6 


non-rad eff(%) 


85.5 ± 1.4 


85.5 ± 1.4 


4f bkg frac (%) 


32.6±1.0 


5.8±1.0 


161 GeV 


Data 


oU4 




200 GeV 


Data 


1636 


1090 


MC Expected 


OQQ -1- A 


97/1 -1- 
Z / 4 It o 


MC Expected 


1717 ±25 


1130±10 


non-rad efF(%) 


83.4 ± 1.0 


80 ± 1 


non-rad eff(%) 


82.6±0.7 


76.9±0.7 


4f bkg frac (%) 


6.2 ± 1.7 


2.1 ± 1.4 


4f bkg frac (%) 


33.6±1.0 


6.0±1.0 


172 GeV 


Data 


280 


zi8 


202 GeV 


Data 


806 


519 


MC Expected 


OSS -1- 7 
Zc5<5 It ( 


9'?9 -1- 
ZoZ It o 


MC Expected 


804±16 


527±5 


non-rad eff(%) 


83 3 ± 1 


79.6 ± 1.0 


non-rad eff(%) 


85.4±0.7 


76.9±0.7 


4f bkg frac (%) 


19.7 ± 2.2 


4.4 ± 1.7 


4f bkg frac (%) 


34.1±1.4 


6.2±1.2 


183 GeV 


Data 


1450 


1 ^77 


205 GeV 


Data 


1687 


1130 


MC Expected 


1478 ± 22 


1083 ±11 


MC Exp(x-t(Hl 


1()93±25 


1089±9 


non-rad eff(%) 


83.0 ± 1.0 


79.0 ± 1.0 


non-rad cff(%) 


82.4±0.7 


76.3±0.7 


4f bkg frac (%) 


27.3 ± 1.2 


5.5 ± 1.2 


4f bkg frac (%) 


34.8±1.0 


6.2±1.0 


189 GeV 


Data 


4448 


3086 


207 GeV 


Data 


2713 


1717 


MC Expected 


4497 ± 37 


3130 ±16 


MC Expected 


2807±32 


1804±12 


non-rad eff(%) 


83.0 ± 0.4 


78.1 ± 0.4 


non-rad eff(%) 


82.7±0.5 


76.5±0.5 


4f bkg frac (%) 


30.0 ± 0.6 


5.5 ± 0.6 


4f bkg frac (%) 


34.7±0.8 


6.2±0.6 



Table 2: Effect of selection cuts (see text) on the number of observed events and the expectation 
from Monte Carlo simulations for the qq signal plus the four-fermion background. The quoted 
errors are purely statistical and do not reflect any experimental or hadronization systematics. The 
efficiency to select non-radiative qq events with — ^fs' < 1 GcV and the expected fraction of 
four-fermion background are also shown. The four-fermion background is negligible below 161 GeV. 



Category 


Error Source 


Standard 


Variation 1 


Variation 2 


(exp.) 


Track-Cluster Matching* 


MT objects 


Tracks + Clusters 




Detector Correction* 


PYTHIA 


HERWIG 




Containment ( | cos | ) * 


< 0.9 


< 0.7 




qqqq (£qqqq) 


< 0.25 


< 0.1 


< 0.4 


qqli^(£qq£l.) 


< 0.5 


< 0.75 


< 0.25 


ISR Algorithm 


Ref. [31] 


Ref. [37] 




Backgrounds (o"bkgd) 


1.0 


+5% 


-5% 


(had.) 


llach'oiiizatioii'' 


PYTIILV 


IIERAMG 




(theo.) 


Fit Range* 


min 6 bins 


+2(1) bins 


—2 bins 


Renorm. Scale Dep. (a;^)* 


1 


0.5 


2.0 



Tabic 3: Summary of systematic variations applied to the datasets. The * indicates which of the 
systematic variations were used for the 91 GeV dataset. 
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91 GeV 


133 GcV 


177 GeV 


197 GeV 




0.1147 


0.1071 


0.0991 


0.0996 


Fit Range (— log^o Vcut) 
/dof 


z.uo— U.oi 

7.98/10 


Z. IO—\J. 10 

7.74/8 


Z. 10— \J. 10 

10.33/8 


z. lO—U. 10 

8.63/8 


Experimental 


±0.0017 


±0.0032 


±0.0019 


±0.0009 


Hadronization 


±0.0026 


±0.0007 


±0.0005 


±0.0005 


Fit Range Variation 

2:^=2.0 
Xfi x^=0.5 


±0.0005 

+0.0030 
-0.0010 


±0.0035 

+0.0030 
-0.0014 


±0.0008 

+0.0023 
-0.0011 


±0.0004 

+0.0020 
-0.0010 


XL=9/4: 


( +0.0040 \ 
I -0.0026/ 


( +0.0033 N 
I -O.OO22J 


( +0.0026\ 
I -0.0018^ 


/ +0.0025^1 
V -0.0019/1 


Theoretical 


±0.0030 


±0.0046 


±0.0024 


±0.0020 


Total Stat. 
Total Syst. 


±0.0004 
±0.0043 


±0.0026 
±0.0056 


±0.0019 
±0.0031 


±0.0008 
±0.0022 


Total Error 


±0.0043 


±0.0062 


±0.0037 


±0.0024 



Table 4: Determination of as and the breakdown of statistical and systematic errors from the fit to 
the Cambridge differential two-jet rate distribution {D2) for all centre-of-mass energy values. The 
quality of the fit is characterized by the chi-square (x^) and the number of degrees of freedom (dof). 
The theoretical error includes contributions from the fit range and variations and excludes the 
XL variation, which is used only for comparison. 





91 GeV 


133 GeV 


177 GeV 


197 GeV 


ds (\/s) 


0.1199 


0.1129 


0.1060 


0.1064 


Fit Range (- logio Vcut) 
xVdof 


1.81-0.68 
4.02/9 


2.75-0.75 
7.02/8 


2.75-0.75 
1.59/8 


2.50-0.75 
9.10/7 


Experimental 


±0.0025 


±0.0026 


±0.0013 


±0.0007 


Hadronization 


±0.0017 


±0.0005 


±0.0005 


±0.0007 


Fit Range Variation 

■^/i x^=0.5 


±0.0005 

+0.0037 
-0.0013 


±0.0020 

+0.0028 
-0.0020 


±0.0035 

+0.0031 
-0.0017 


±0.0015 

+0.0030 
-0.0014 


Xi=9/4 

■^J- XL=4:/9 


/ +0.0044 \ 
I -O.OO27J 


/ +0.0028\ 
I -0.0026J 


/ +0.0031 \ 
I -O.OO22J 


/ +0.0032-, 
>^ -0.0020; 


Theoretical 


±0.0037 


±0.0034 


±0.0047 


±0.0034 


Total Stat. 
Total Syst. 


±0.0004 
±0.0048 


±0.0026 
±0.0043 


±0.0019 
±0.0049 


±0.0009 
±0.0035 


Total Error 


±0.0048 


±0.0050 


±0.0052 


±0.0037 



Table 5: Determination of as and the breakdown of statistical and systematic errors from the fit 
to the Durham differential two-jet rate distribution (D^) for all centre-of-mass energy values. The 
quality of the fit is characterized by the chi-square (x^) and the number of degrees of freedom (dof). 
The theoretical error includes contributions from the fit range and Xf^ variations and excludes the 
XL variation, which is used only for comparison. 
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91 GeV 


133 GeV 


177 GeV 


197 GeV 




0.1254 


0.1158 


0.1064 


0.1066 


Fit Range (- logio ^cut) 


2.56-0.56 


2.75-0.50 


2.75-0.50 


2.75-0.50 




20.70/16 


7.38/9 


2.57/9 


5.16/9 


Experimental 


±0.0013 


±0.0040 


±0.0030 


±0.0015 


Hadronization 


±0.0033 


±0.0034 


±0.0018 


±0.0017 


Fit Range Variation 


±0.0008 


±0.0042 


±0.0020 


±0.0010 


:r,,=2.0 
•'■/' ,c„=0.r) 


+0.0036 
-0.0001 


+0.0021 
-0.0001 


+0,0016 
-0.0001 


+0,0016 
-0.0001 


Theoretical 


±0.0037 


±0.0047 


±0.0026 


±0.0019 


Total Stat. 


±0.0005 


±0.0033 


±0.0023 


±0.0009 


Total Syst. 


±0.0051 


±0.0070 


±0.0042 


±0.0030 


Total Error 


±0.0051 


±0.0077 


±0.0048 


±0.0031 



Table 6: Determination of ctg and the breakdown of statistical and systematic errors from the fit to 
the Cambridge average jet rate distribution {{N)^) for all centre-of-mass energy values. The quality 
of the fit is characterized by the chi-square (x^) and the number of degrees of freedom (dof). 





91 GeV 


133 GeV 


177 GeV 


197 GeV 




0.1272 


0.1193 


0.1103 


0.1106 


Fit Range (- logio Vcut) 


2.56-0.56 


2.75-0.50 


2.75-0.50 


2.75-0.50 


xVdof 


7.90/16 


4.73/9 


0.81/9 


4.14/9 


Experimental 


±0.0006 


±0.0039 


±0.0026 


±0.0013 


Hadronization 


±0.0039 


±0.0037 


±0.0009 


±0.0007 


Fit Range Variation 


±0.0003 


±0.0029 


±0.0003 


±0.0008 


™ ^^=2-0 
■''M a;^=0.5 


+0.0034 
-0.0013 


+0.0028 
-0.0010 


+0.0023 
-0.0008 


+0.0023 
-0.0008 


Theoretical 


±0.0034 


±0.0040 


±0.0023 


±0.0024 


Total Stat. 


±0.0004 


±0.0030 


±0.0021 


±0.0008 


Total Syst. 


±0.0052 


±0.0067 


±0.0035 


±0.0028 


Total Error 


±0.0052 


±0.0073 


±0.0041 


±0.0029 



Table 7: Determination of and the breakdown of statistical and systematic errors from the fit to 
the Durham average jet rate distribution {{N)^) for all centre-of-mass energy values. The quality 
of the fit is characterized by the chi-square (x^) and the number of degrees of freedom (dof). 
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91 GeV 


Dk 


^2 






133 GeV 












100 


85 


53 


46 


^2 


100 


83 


70 


64 




85 


100 


56 


47 


^2 


83 


100 


64 


63 




53 


56 


100 


75 




70 


64 


100 


81 




46 


47 


75 


100 




64 


63 


81 


100 




177 GeV 


Di' 




(NY 




197 G(V' 




Di' 


(Nf 




^2 


100 


84 


73 


66 


^2 


100 


90 


77 


68 


^2 


84 


100 


71 


72 


^2 


90 


100 


72 


73 




73 


71 


100 


82 




77 


72 


100 


82 




66 


72 


82 


100 




68 


73 


82 


100 



Table 8: Statistical correlation matrix determined for the combination of as measurements from the 
four different observables for each of the four different datasets. Each matrix element is presented 
as a percentage (%). 



r~s (GeV) 




Value for ^ 


/i 




Value at Mz 




Cstat. 




CThadr. 


•^theory 


Us 


Cstat. 




CThadr. 


•^theory 


Mz 


0.1213 


0.0004 


0.0013 


0.0029 


0.0034 


0.1213 


0.0004 


0.0013 


0.0029 


0.0034 


133 


0.1126 


0.0025 


0.0028 


0.0007 


0.0039 


0.1191 


0.0028 


0.0031 


0.0008 


0.0044 


177 


0.1039 


0.0018 


0.0018 


0.0001 


0.0028 


0.1140 


0.0021 


0.0021 


0.0001 


0.0033 


197 


0.1046 


0.0008 


0.0009 


0.0002 


0.0023 


0.1163 


0.0010 


0.0012 


0.0003 


0.0029 



Table 9: Values of ag determined from the weighted average of the individual CKg results at each 

centre-of-mass energy, along with statistical and systematic errors. The value for the dataset when 
run back to the Z pole are given on the right half of the table with full statistical and systematic 
errors. 
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Figure 1: The y/s — Vs' distribution for the 189 GeV dataset. The vertical hue indicates where the 
cut was apphed. The expected non-radiative qq signal — V^'^true ^ GeV) was determined from 
Monte Carlo and normalized to the luminosity of the measured sample. Vertical error bars indicate 
the size of statistical errors, while horizontal error bars correspond to the bin width. 
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OPAL (189 GeV) 




qqqq 



OPAL (189 GeV) 



Total Expected 
qq Signal 

4f(qqll)_ 

4f(qqqq) 

Data 
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0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



0.9 



L 



qqlv 

Figure 2: Distributions of the four-quark likelihood >Cqqqq (top) and semi-leptonic likelihood >Cqq^^ 
(bottom) for the 189 GeV dataset. The expected contribution of the four-quark non-QCD back- 
ground to the total sample is shown by the dotted line, and the expected semi-leptonic four-fermion 
background is represented by the dashed line. The vertical lines indicate the values of the pa- 
rameters where the cuts were applied. The expected qq signal was determined from Monte Carlo 
and normalized to the luminosity of the measured sample. Vertical error bars indicate the size of 
statistical errors, while the horizontal bars correspond to the bin width. 
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Figure 3: The hadron level n-jet rates for the R and e variants of the Cone algorithm for the data 
with x/s = 91 GeV (top) and y/s = 133 GeV (bottom). In all plots the JETSET/PYTHIA and 
HERWIG Monte Carlo expectations arc represented by the curves. Outer error bars indicate total 
errors while the inner bars indicate statistical errors. 
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Figure 4: The hadron level n-jet rates for the R and e variants of the Cone algorithm for the data 
with ^ = 179 GcV (top) and = 198 GcV (bottom). In all plots the PYTHIA and HERWIG 
Monte Carlo expectations are represented by the curves. Outer error bars indicate total errors while 
the inner bars indicate statistical errors. 
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^cut ^cut 



Figure 5: The hadron level n-jet rates for the JADE algorithm for the data with y/s = 91 GeV (top 
left), ^/s = 133 GeV (top right), ^/s = 179 GeV (bottom left) and ^/s = 198 GcV (bottom right). 
In all plots the JETSET/PYTHIA and HERWIG Monte Carlo expectations are represented by the 
curves. Outer error bars indicate total errors while the inner bars indicate statistical errors. 
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Figure 6: The hadron level n-jet rates for the Durham and Cambridge algorithms for the data with 
^/^ = 9l GeV (top) and ^ = 133 GcV (bottom). In all plots the JETSET/PYTHIA and HERWIG 
Monte Carlo expectations are represented by the curves. Outer error bars indicate total errors while 
the inner bars indicate statistical errors. 
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Figure 7: The hadron level n-jet rates for the Durham and Cambridge algorithms for the data with 
^ = 179 GeV (top) and ^ = 198 GeV (bottom). In all plots the PYTHIA and HERWIG Monte 
Carlo expectations are represented by the curves. Outer error bars indicate total errors while the 
inner bars indicate statistical errors. 
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Figure 8: The differential 2-, 3- and 4-jet rates at the hadron level as a function of ycut for the 
Cambridge and Durham algorithms for data with = 91 GeV (top) and for the LEP1.5 com- 
bined dataset at y/s = 133 GeV (bottom). PYTHIA and HERWIG Monte Carlo expectations are 
represented by the curves. The differential 3- and 4-jet rates have a negative slope in the region of 
large ycut and therefore the negative values of D3 and D4 are plotted (the positive values are not 
seen on the curves). Note that D3 and D4 are scaled down by one and two orders of magnitude, 
respectively, for clarity. Error bars indicate total (statistical -|- systematic) errors. 
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Figure 9: The differential 2-, 3- and 4-jet rates at the hadron level as a function of ycut for the 
Cambridge and Durham algorithms for two combined LEP2 datasets with =177 GeV (top) 
and y/s = 197 GeV (bottom). PYTHIA and HERWIG Monte Carlo expectations are represented 
by the curves. The differential 3- and 4-jet rates have a negative slope in the region of large 2/cut 
and therefore the negative values of Ds and D4 are plotted (the positive values are not seen on the 
curves). Note that Ds and D4 are scaled down by one and two orders of magnitude, respectively, 
for clarity. Error bars indicate total (statistical + systematic) errors. 
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Figure 10: The average jet rates at the hadron level as a function of ycut for the Cambridge algorithm 
for the all centre-of-mass energies. PYTHIA and HERWIG Monte Carlo expectations are represented 
by the curves. Error bars indicate total (statistical + systematic) errors. 
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Figure 11: The average jet rates at the hadron level as a function of ycut for the Durham algorithm for 
the all centre-of-mass energies. PYTHIA and HERWIG Monte Carlo expectations are represented 
by the curves. Error bars indicate total (statistical + systematic) errors. 
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Figure 12: Fits of the parton level D2 distribution using the Cambridge algorithm as a function of 

ycut to the In R prediction for the 91 GcV (top left), 133 GeV (top right), 179 GeV (bottom left) 
and 198 GeV (bottom right) datascts. The inset plots show the pull of each point from the line 
of best fit vs log^oJ/cut- The two vertical lines delineate the fit range. Vertical error bars represent 
statistical errors only. 
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Figure 13: Fits of the parton level D2 distribution using the Durham algorithm as a function of 

Vent to the In R prediction for the 91 GeV (top left), 133 GeV (top right), 179 GeV (bottom left) 
and 198 GeV (bottom right) datascts. The inset plots show the pull of each point from the line of 
best fit vs log^oJ/cut- The two vertical lines delineate the fit range. Vertical Error bars represent 
statistical errors only. 
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Figure 14: Fits of the parton level (N) distribution using the Cambridge algorithm as a function of 
ycut to the In R prediction for the 91 GeV (top left), 133 GeV (top right), 179 GeV (bottom left) 
and 198 GeV (bottom right) datasets. Vertical error bars represent statistical errors only. Note: 
pulls are not shown as there are very large bin-to-bin correlations. 
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Figure 15: Fits of the parton level {N) distribution using the Durham algorithm as a function of 
ycut to the In R prediction for the 91 GeV (top left), 130 GeV (top right), 189 GeV (bottom left) 
and 207 GeV (bottom right) datasets. Vertical error bars represent statistical errors only. Note: 
pulls are not shown as there are very large bin-to-bin correlations. 
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Figure 16: Comparison of combined ctg values with those determined separately from the D2 , D)^, 
(N)'^' and (N)^ distributions for each centre-of-mass energy dataset. Outer error bars indicate the 
size of the total errors while inner bars indicate the size of the statistical errors. 
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Figure 17: Comparison of the ag values after having been run back to the Z pole for each of the 
datasets. The point labelled by 'ALL' represents the value determined from the weighted mean 
of the four combined as{Mz) determinations. The shaded band corresponds to the one standard 
deviation range of the world average value of as(-^z) [41]- The inner error bars represent the 
statistical errors while the outer error bars represent the total error. 
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Figure 18: The running of ag as a function of centre-of-mass energy. The points correspond to the 
Z-cahbration dataset, the combined LEP1.5 data and the two LEP2 regions 161-183 and 189-209, 
respectively. The soUd hue corresponds to the expectation based on the weighted mean of the four 
combined as{Mz) determinations, and the outer dashed hne to the total (statistical + systematic) 
uncertainty. 
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